clc;close;clear;
X = [8, 4, 5; 4, 8, 6; 5, 6, 9];
% tic
% [L, p] = Cholesky(X);
% toc
tic
b = chol(X,'lower');
toc
tic
[L, p] = Cholesky(X);
toc
disp(L)
disp(b)